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The kaon-antikaon system is studied in both the atomic and the strongly interacting sector. We 
discuss the influence of the structures of the /o(980) and the ao(980) mesons on the lifetime of 
kaonium. The strong interactions are generated by vector meson exchange within the framework 
of the standard SU(3)v x SU(3)a invariant effective Lagrangian. In the atomic sector, the energy 
levels and decay widths of kaonium are determined by an eigenvalue equation of the Kudryavtsev- 
Popov type, with the strong interaction effects entering through the complex scattering length for 
KK scattering and annihilation. The presence of two scalar mesons, /o(980) and ao(980), leads to 
a ground state energy for the kaonium atom that is shifted above the point Coulomb value by a few 
hundred eV. The effect on the lifetime for the kaonium decay into two pions is much more dramatic. 
This lifetime is reduced by two orders of magnitude from 1.2 x 10 -16 sec for annihilation in a pure 
Coulomb field down to 3.2 x 10~ 18 sec when the strong interactions are included. The analysis of 
the two photon decay width of the fo (980) suggests a generalization of the molecular picture which 
reduces the lifetime of kaonium still further to 1.1 x 10 -18 sec. 
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I. INTRODUCTION 

There has been substantial experimental progress in the field of meson spectroscopy during the last decade [lj. In 
the energy region up to 2 GeV, more scalar-isoscalar mesons have been established than can be accounted for by a 
quark-antiquark structure 0. The structure of the scalar meson with lowest mass, the /o(980), has been controversial 
for many years. The / (980) might be a qq state 0,0, a q 2 q 2 state 0, or a KK molecule 0. 

The radiative decay of the (p-meson provides a particularly strong argument for a q 2 q 2 interpretation of both 
the /o(980) and the ao(980) mesons, as was first pointed out by N.N. Achasov and V.N. Ivanchenko,5J. Both the 
recent Novosibirsk dataQ] and the KLOE data0 can be reproduced by a calculation which generates those mesons 
dynamically, however 0. Since in Ref.0 Oiler had to introduce a 4> — > ^K°K contact interaction, the issue remains 
controversial ^3 • 

The production of two neutral pions in ultrarelativistic pion proton reactions shows a strong dependence of the 
S-wave amplitude on the momentum transferred between the proton and the neutron for invariant two-pion masses 
in the vicinity of 1 GeV0,^2- This fact has been interpreted as evidence for a hard component of the /o(980), see 
e.g. Ref. 01 . A recent calculation based on a model which allows for a dynamical generation of the /o(980) provides 
a good description of the data^4|. On the other hand, the model employed in Ref.[L"H also includes a bare scalar 
resonance. Given this situation, we feel that a simplified calculation might be helpful. 

In this paper, we develop an analytical model which generates the /o(980) meson as a bound KK structure. This 
model makes specific predictions for the structure of the exotic atom Kaonium. In a second step, we work out the 
predictions for Kaonium based on the meson-exchange model of Ref. . 

The molecular interpretation is consistent with the small binding energy ~ 10-20 MeV of the KK system relative to 
its reduced mass Mk ~ 496 MeV. This sugg ests a non-relativisic effective field theory approach 15] that has also been 
used rece ntly to study both pionium[l6l. Il7| as well as 0(a 6 ) QED recoil and radiative corrections to the positronium 
spectrum^. With this in mind we use the standard SU(3)v x SU(3)a Lagrangian to describe the dynamics 
of the KK interaction pfl l2l| and decay via the exchange of p, u, <j>, K*.... vector mesons. The coupling constants 
appearing in the Lagrangian are related via SU(3) symmetry to the p-n-n coupling constant gpn-n, which in turn can 
be obtained from the KSRF relation 22J. Thus once the decision on the form of Lagrangian has been taken, only 
physical meson masses and known coupling constants enter into the calculations. 

We work in the non-relativistic limit in which case two important simplifications occur: (i) only r-channel exchange 
diagrams survive, and (ii) the resulting one-meson exchange potentials become local. This in turn means that one can 
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reduce the Bcthe Salpeter (BS) integral equation for the bound states of the interacting KK system to a local two 
body Schroedinger wave equation that offers a significant simplification over working with integral eauations |20ll23| . 

After a brief recall of the derivation of the wave equation for non-relativistic local potentials from the BS equation 
in Section II, the calculation of one meson exchange potentials involving both direct p, w, </> transfer between K and K 
as well as KK — > KK scattering via two-pion intermediate states involving K* exchanges is carried out in Section III. 
We refer to these potentials collectively as one-boson exchange (OBE) potentials. The last contribution is essential to 
describe the KK — > 2tt decay channel. Then we make use of the fact that one can describe the low-energy properties 
of KK very adequately in the effective range approximation, to replace the OBE potentials by phase-equivalent 
potentials of the Bargmann type[24| that give rise to the same scattering length and effective range. These potentials 
offer the unique advantage of having known analytic solutions so that the associated Jost functions can be constructed 
explicitly. 

A knowledge of the Jost functions in turn determines both the scattering and bound state properties of the KK 
system without futher approximation 25j. In Section IV we carry out this program and compute both the mass and 
decay width of the kaonic molecule from the relevant Jost function that includes annihilation contributions. The 
resulting complex total energy for this system is (981 — 25z) MeV that is in reasonable agreement with the recent 
experimental data from Fermilab|26| that give [(975 ± 3) - (22 ± 2)i] MeV. 

We also give computed elastic and reaction cross sections for KK scattering, the irir inelasticities, and the cross 
section for the inverse process, tt + it~ — ► KK, for which data exists, using detailed balance. Those calculations 
disagree with the measured cross section, particularly near threshold |27t 1281 129| . 

On the other hand the similarity in the pole position in both the recent Fermilab measurements as well as the 
calculated position of this pole is striking. We also show that the molecular picture is totally inadequate for describing 
the 2 photon decay width of /o(980), which comes out to be an order of magnitude larger than experiment. Taken 
together with the overall underestimate of the ir + ir~ — > KK data, this suggests that the /o(980) ground state cannot 
be a pure molecular state (see also [13 )• This point is taken further at the end of Section ITvl 

We take up the discussion of kaonium in Section V. Since kaonium is a mixture of isoscalar and isovector states due 
to isospin-breaking introduced by the Coulomb potential, we appeal to the no-internal-mixing approximation that 
was introduced in connection with isobaric analog states in nuclei. In this approximation one joins linear combinations 
of the I — and 1=1 isospin amplitudes in the external region where only the Coulomb field is important, onto 
(in this case) known solutions of good isospin of the wave equation in the internal region where the strong field is 
dominant. This allows one to construct the Jost function for kaonium without further approximation. Only momenta 
of O(Mko), where a « 1/137, are of interest for kaonium. This fact allows one to recast the condition that its bound 
states are given by the zeros of the Jost function in the lower half of the complex momentum plane as an eigenvalue 
equation of the Kudrvavtsev-Popov tvpe|3lj . in which the strong field effects only enter through the K + K~ scattering 
length. The solutions of this equation show that the kaonium levels are shifted upwards from their pure Coulomb 
values and acquire lifetimes between 10 _18 -10 -16 sec in the presence of the KK molecular ground state. 

We also show explicitly that the strong field introduces an additional node into the eigenfunctions of kaonium 
that lies at the K + K~ scattering length. The kaonium levels are therefore to be viewed as excited states of the 
kaonic molecule in the combined Coulomb and strong fields of the system. This in turn has the effect of enhancing 
the K + K~ — > 2tt annihilation widths by two orders of magnitude over what they would have been under Coulomb 
binding of the charged kaon pair alone. 



II. FORMALISM 



In order to obtain a non-relativistic equation for describing possible bound states of a KK meson pair resulting 
from vector meson exchange between them we briefly recall the derivation given by Landau and Lifshitz 32] that starts 
out with the four-point vertex ir[p3,p4;pi,p2] that enters the Bethe-Salpeter (BS) equation^]. Thepi,p2 and^3,p4 
are incoming and outgoing meson four momenta respectively. Then the homogeneous equation for this vertex that 
determines the bound state poles of the BS equation is 

^^zf [p 3 ,p 4 ;g,P3 +P4 - q](— i)D(q)(-i)D(p 3 + p 4 - q)iT[q,p 3 +p 4 - q\P\P2\ (1) 

All momentum labels formally flow from right to left, and the sum of each pair on either side of the semi-colon equals 
the total four momentum P of the incoming pair which is conserved throughout the diagrammatic equation. The 
hatted vertex is the irreducible piece that generates T by iteration. The D's are meson propagators. For a free meson 
of four momentum q and mass M one has 

_ = q 2 - M 2 + iO ~ ^2m\ -M -q 2 /2M + iO ® 
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if in addition we move to the non-relativistic limit by formally suppressing propagation "backwards" in time, i.e. 
by omitting the anti-particle pole in the upper half of the complex q plane. We note[32j that the momenta P\,P2 
are simply labels in that are not determined by the equation at all. So one can simply drop them. The other 
simplification is to observe that the combination 

x(P3,Pa) = D(p 3 )T[p 3 ,p4]D(p 4 ) (3) 

appears under the integral sign. Thus one can equally well recast Eq. JJJ as an integral equation for xiPZiPt)- This 
is most usefully written down in the center of mass (CM) system for which [^3,^4] = ±p + P/2 with p — [po, p] and 
P = [Po, 0], where Po is the total energy in the CM frame. Then the equation determining the bound states reads [3^ 

i X {p, P) = D{p + \p)D{-p + lp) j -0jf \p + ip, -p + \P\q + \P, -q + \p]x{q, P) 

(4) 

where x(Pa,Pi) = x(p + P/2, -p + P/2) = X (p, P)- 

To make further progress towards a non-relativistic equation, the vertex T should not depend on the time compo- 
nents of the relative outgoing and incoming four momenta p = [p , p] and q — [q 0l q], i.e. 

f[ P + ip, -p +\p;q + \p, -q + \p] « f [ P , q , p ] (5) 

Should this be the case one can then integrate out these time components from \ to obtain what is effectively the 3D 
fourier transform of the wave function for relative motion, 

iKp)= r ^r^ p ) ( 6 ) 

Carrying out the integrals over the time components of the relative four momenta one finally arrives at the desired 
equation, 

,P 2 , OJf u\if \ f d 3 q f[p,q,P ] , . 

(- + 2M PoMp) - j {2jr)3 ^(q) = (7) 

that makes use of the non-relativistic approximation 

d £ D(P + '-mm-r + \p») ( Wl + \ u _ Po ) (8) 

as given by Eq. . One recognises Eq. J7J as the Schroedinger wave equation in momentum space for relative motion 
in a potential 

v(P, q) = -%5! (9) 

and binding energy e = 2M — P n > 0. Should T[p, q, P ] only depend on the difference p — q, the corresponding 
potential will be local in coordinate space. Since T is a relativistic proper vertex that in lowest order gives the T 
matrix, the mass factor 1/4M 2 that converts f into a non-relativistic potential is the same factor that relates the 
relativistic T matrix to its non-relativistic counterpart. 



III. ONE BOSON EXCHANGE POTENTIALS 



In the following we use Eq. 10 to investigate the possible binding of the KK system via one-boson exchange 
potentials by constructing T from the relevant pieces of the standard SU(3)v X SU(3)a invariant Lagrangian, the 
derivation and properties of which are described in detail in Rcf. 19] . For KK scattering the relevant interaction 
Lagrangians for our purposes are 

Ckk p = gKKplK^rd^K}^ + cc 
Ckklo = gKKujiK^d^K}^ + cc 

l~KK(t> = gKK4> [K^d^K]^ + cc (10) 
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which, together with 

£>ttKK* = 9-nKK* [0**7f • K^tK; - 7? • d^K^rK;} + cc (11) 

generate interaction potentials for KK — > if if scattering via p, w, vector meson exchange, as well as for if if — > 7T7t 
annihilation via if* (892) strange meson exchange. The K's are all isospin doublets and cc stands for the additional 
charge conjugation contribution with if — > AT, if* — > if*. 

The coupling constants in these expressions are all fixed in terms of the pirn coupling constant gp-mx by 517(3) 
symmetry relations|l9|. 

1 

gKKp — 9kk p — 2^P W7r 

1 

gKKtu — 9KKu> — ~ 2 9 pirn 
1 

9KK4, = 9kR4> - ~~^ g P* 7T 

1 

9-kKK* = 9-kRk* - ~2 g P™ ( 12 ) 

On the other hand the p7T7r coupling is determined by the KSRF relation as g p7r7T » M p /(V2f 7I ) « 6 in terms of the 
p meson mass M p and the pion weak decay constant / ff « 93MeV. In this sense, then, there are no free parameters in 
the calculation of the exchange potentials. They only contain physical meson masses and known coupling constants. 

In the non-relativistic limit, pi + pj [Mi + Mj , 0] , p; — pj ss [0, pt — Pj] only the t-channcl scattering diagrams 
are relevant for determining T. In this limit these amplitudes all have a common Yukawa like form 

1 1 1 1„, 2 ^ (Mi + M 3 )(M 2 + MA , , 

nv + - 2 P, p +- 2 P; q + - 2 P, , + - 2 P] « ?Cn > 41 (13) 

in momentum space, where k = p — q is the t channel 3-momentum transfer, and Mi are the masses in the indicated 
entrance and exit channels; M is the mass of the exchanged boson and g the coupling constant. The isospin and 
boson identity factors Cj and n s in Eq. I|13|l are given in|2(j. They are Co = [3, 1, 1] and C\ = [—1, 1, 1] and n s = 1 
for p, u>, <f> exchange in the if K — > if K isoscalar and isovector channels. The corresponding values for if* exchange 
in the if K — > tttt isoscalar and isovector channels are Co = — -y/6, C\ = —2, and n s = 1/V2- 



A. ifif exchange potentials 

Since T is only a function of the 3-momentum transfer, the resulting potentials in Eq. 10 are attractive local Yukawa 
potentials in coordinate space, and the determination of their possible bound states is numerically straightforward. 
While not necessary for convergence in the non-relativistic case, we, however, also include at each vertex in Eq. 113|) 
the form factor containing an arbitrary cutoff A 



F*(k) = 



2A 2 - M 2 \2 



2A 2 + k 2 



(14) 



that was used in |2flj, [23| to obtain convergence in scattering calculations based on the relativistic T vertices. ^From 
Eq. 10 the coordinate space non-relativistic if K — > if K potential associated with the exchange of boson M is then 



V(r) = -g 2 C I U{M, r) 

after division by 4M| f where Mk is the kaon mass, and U(M,r) is the fourier transform 



(15) 



U(M, r) = 



d 3 k [F'(k)] : 



ik-r 



(2tt) 3 M 2 + k 2 



4tt I 



-Mr 



-V2Ar 



1 , M 2 Ai 4 w M 2 W - s 



— (3 77?)(1 r^7 (V2Ar) 2 + — 1 — 3 v^Ar) 3 ^ 

16 v 2A 2A 2A 2 ^ v 7 48 v 2A 2 ^ v '\\ 



(16) 
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The resulting K K potential due to p, u>, <fi exchange are attractive in both isospin channels 7 = and 1, 

V I=Q (r) = -9%„{\TJ{M p ,r) + \u{M u ,r) + \u{M^r)) (17) 

and 

7/=i(r) = -g% n { - \u{M p ,r) + \u{M^,r) + \u{M^r)) (18) 

but too weakly so in the latter channel due to an almost complete cancellation between p and lo exchange to support a 
bound state. In deriving these expressions we have exploited the SU(3) relations l|12f) to express all coupling constants 
in terms of <7 p7r7r . 

These potentials are plotted in Fig. [J for a strong coupling constant a s = g 2 plT7T /AiT = 2.9 and cutoff A = 4 GeV. 
In principle A can be different for each exchanged boson[2(|. One notes that the Vj are now finite at the origin, in 
contrast to the sum of pure Yukawa potentials to which they reduce as A — > oo. 

The coordinate space version of the wave equation Eq. Q) that describes the relative s-state motion in these 
potentials becomes the usual Schroedinger equation 

[V 2 + fc 2 - VKr)hKr) = (19) 

with total CM energy Pq — 2Mk + k 2 /Mk- We ignore the charged to neutral kaon mass difference and work with an 
average kaon mass Mk ■ By numerically constructing the k — > scattering solutions one easily obtains the scattering 
lengths and effective ranges that enter into the effective range expansion of the s-wave phase shifts S(k) 

k cot S(k) = — + - ri k 2 + ■■ ■ (20) 
aj 2 

in the two channels as 

a = 5.835M K \ r = LINM' 1 (isoscalar) 
ai = -0.401M^\ n = 3.702Af K 1 (isovector) (21) 

Thus only the isoscalar potential supports a bound state. Its binding energy is calculated to be e = 18.63 MeV, 
placing the total mass for the bound KK system at Pq = 2Mk~£ = 973.4 MeV. Fig. [21 illustrates the sensitivity of e, 
ao and ro on the choice of cutoff for this channel. For A > 4 GeV the values of do and ro are only weakly dependent 
on A. 

Even though the isovect or p otential does not support a bound state the ao(980) can be fully explained by this 
model as a threshold effect |2l|. Effects of the ao(980) thus are included in the calculation though it should be noticed 
that decays to the 7r^-channel are not considered in the analytical model. 



B. Shape— independent approach 

The relatively small (on a hadronic scale) binding energy e/M^ ~ 0.04 in the isoscalar channel (and no binding 
at all in the isovector channel) indicates that an effective range approach should be applicable. This is confirmed 
by recalculating the binding energy from the effective range formula pEj e « k 2 /Mk — 18.60 MeV where k = —a = 
0.194A/R-, see Eq. (|27ll below. Consequently it is entirely sufficient to characterize the KK interaction at low energies 
in terms of a scattering length and an effective range: the actual shape of the potential is immaterial. 

We make use of this feature in the following to replace the strong potentials Vi in Eqns.ll7|) and (11811 by analytic 
potentials of the Bargmann type given below that are phase equivalent to them. Bargmann|24| (see also 33] for an in 
depth overview) has shown how to construct families of potentials that give rise to a prescribed Jost function f(k). 
If one chooses 

/(*) - (22) 



then the potential that leads to this Jost function can be determined. It is 50] 24, 3 



1 8& 2 r e br 



M K b 2 - a 2 



b- 



Vi(r) = -TT-7T— » 7— + 7— (23) 
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For this potential the effective range expansion is exact: 



k cot S(k) 



ba 



b — a b — a 



This allows one to identify the scattering length and effective range as 

b- a 2 



a I 



ba ' 



ri 



b-. 



(24) 



(25) 



where the actual values of the parameters a and b will depend on the isospin channel. 

Fixing a and b is obviously equivalent to prescribing the scattering length and effective range. If a is negative, 
a = —k, say, then f(k) has a zero at k — —in in the lower half of the complex k plane, and there is a bound state |33j 
of binding energy e s = k 2 /Mk- Both the scattering and bound state wave functions of the Bargmann potential 
in Eq. 12- il) are known explicitly. Further details will be found in the Appendix. The availability of such analytic 
solutions will be important for including the effects of the K K — > tttt decay channel non-perturbatively as discussed 
below. 

Now invert Eq. I|25[) to construct phase equivalent potentials that have the same scattering length and effective 
range as the original exchange potentials in Eqs. (fT7| and fTTfy . The required values of a and b are 



1 r 

ri 

b= 1 
ri 



1 - y/1 - 2r//a/ 



1 + y/l- 27-// 



a i 



or 



a = -0.1936M K 
a = +0.9219M K 



b = 1.491M K 
b= 1.462M K 



(isoscalar) 
(isovector) 



(26) 



(27) 



for the two channels in question. These parameters are summarized in Table[l]for easy reference. The resulting set of 
phase equivalent potentials are included in Fig. ^ 



C. KK annihilation potential 
From Eq. Ijl3|) one reads off the irreducible vertex for KK — > tttt via K* exchange, as 



2 (Mtt + M K 



,2 



T(KK - tttt) a -Cm s gi KK , v * (28) 



Ml, 

where n s Ci = —^/3 or — \[2 in the isoscalar or isovector channel. We have neglected the 3-momentum transfer relative 
to the large if* (892) mass. As shown in Fig. [3 this amounts to introducing point vertices that can contribute to KK 
scattering via an intermediate pion loop. The contribution from this interaction is complex due to the possibility of 
on-shell decays KK — > tttt in the intermediate state, and the resulting KK bound state acquires a width. The details 
are as follows: calling the contribution from Fig. |3«r^ 1 - ) one finds 

= it*(KK _ ^ / (03 12/M7r + 2 MT -P o 

= if\KK^™)^ (29) 

in view of Eq. iJSJ. As summarized briefly in Eqs. (|A.23() and (|A.24|I of the Appendix the three dimensional integral 
J(Pq) can be dimensionally regulated in d dimensions without difficulty to obtain a finite result |16| for d — > 3. The 
associated annihilation potential V K ^ is then obtained from Eq. ^ as an attractive contact potential in coordinate 
space, 

VKR{r) = - w^pir J(Fo) ' 3(r) = ~ c ^ 3(r) (30) 
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with 



3 (91kk») 2 fM*+M 



16 M 2 M 2 K 



The value of this factor in the isovector channel is c 2 — |cg. The function 



J{Po) 



.Mtt_ r 

' 47T 



M 7r {P Q -2AU] 



1/2 



(31) 



(32) 



is a complex function of Pq that is positive and pure imaginary along the upper lip of the branch cut extending from 
2M„ to oo along the real Pq axis, see Eq. (|A.24ll . 



D. Jost function for a finite plus delta function potential at the origin 

We investigate the influence of the annihilation potential on the properties of the KK system by constructing the 
revised Jost function for the sum Vj + V K g. Since Vkk i s a contact interaction, and the solutions in the potential Vi 
are known, this can be done without approximation. The new scattering problem to be solved reads 

[V 2 + k 2 - M K Vi{r) + M K 48 3 (r)]ip(r) = (33) 

Consider isospin 1 = 0. The presence of the delta function potential obliges one to give up the boundary condition 
ip(k,r) ~ r — > on the scattering wave function tp(k,r) = rip(k,r) at the origin in favor of prescribing j^lj its 
logarithmic derivative y'/y there |33|. By integrating Eq. i.Y.'A) over the volume of an infinitesimal sphere centered 
at the origin, one finds that the logarithmic derivative at the origin is fixed by the strength of the delta potential 
according to 

The point is now that for r > both the boundary condition on (p(k, r) and the closed form of irregular solution of 
Eq. (|33|l are known explicitly: they are given by Eqs. I|34(l and (|A.4(l respectively. By constructing the Wronskian of 
this pair in the limit r — > one obtains the revised Jost function as 

f(k) = W[f,v]=\(( + ik) } ^ + i b -^j-] (35) 

after setting lim ip(k,r) — > 1 for convenience. This choice is unimportant since neither the root of f(k) nor the 
scattering phase shift depends on it. The value of the logarithmic derivative C is k— dependent through P : from 
Eqs. (PU and (Hlj! one finds 

c«-^«-(35Tk)V(ft-^)]-^ 

P = 2M K + ^— (36) 

M-K 

after taking the SU(3) value g^pcx* = — Qp^-nl^ for the coupling constant. 



IV. NUMERICAL RESULTS 



The single equation 135|) contains all the necessary information on the scattering and bound states of the two-body 
K K system. Moreover, this information can be extracted without further approximation. We use this form of f(k) to 
identify the revised scattering length and effective range as well as the bound state energy in the presence of KK — > 2ir 
annihilation. 
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1. Effective range expansion 



The KK —f 7T7T annihilation channel renders the scattering matrix S{k) = f(k)/f(—k) non-unitary, since now 
f*(k) f(—k). This can be seen directly from Eq. I|35|) after noting that the logarithmic derivative C is even 
in k and purely negative imaginary for real k. The scattering phase shift is now complex. In order to identify the 
scattering length and effective range in the presence of absorption, one can either maintain the expansion (|20() , thereby 
introducing complex effective range parameters |30l l3l| . or define [34[ these in terms of the real part S(k) of the total 
phase shift (not the phase of S(k)), 

S(k) = f(k)/f(-k) = V e 2 ^ k \ < r? < 1 (37) 

In the latter case 

t *X(* . J(k)\f(-k)\ + f(-k)\f(k)\ 

k cot 6(k) = lkm - = ^ fm{ _ k)] _ f{ _ k)][m (38) 

that shows explicitly that k cot S(k) is an even function of k. Set £(fc) = —i£,{k) in Eq. I|35[l for the Jost function, and 
then expand the right hand side to 0(k 2 ). After some calculation one establishes the revised effective range expansion 

k cot 5{k) = -^r + lr' k 2 + ■■■ (39) 
a 2 

with the exact expressions 

16; a 2 e + (b 2 -a 2 ) 2 \ 



i' b-aiaC 2 - 2ab 2 + a 3 - b 3 . 



(40) 



and 



1 , _ a 2 (b 2 - a 2 ) 4 + a(b + a)[a 3 (b - a)(5b 2 - 3a 2 ) + o 4 (a 2 - 2ab + 2fo 2 )]g 2 
2 r ° ~ (b - a)[a 2 e + {b 2 - a 2 ) 2 ] [a 3 - b 3 - lab 2 + a£ 2 ] 2 + 

a(b + a)[a(b - a)(b 2 - 3a 2 ) - o 4 ]£ 4 + a 4 £ 6 
(b - a)[a 2 £ 2 + (b 2 - a 2 ) 2 ] [a 3 - b 3 - 2ab 2 + a£ 2 ] 2 

4 [(a 2 + ab - b 2 ){a 3 ab 2 b 3 ) + a 2 £ 2 (b + a)] 
C? (o-a)[a 2 £ 2 + (o 2 -a 2 ) 2 ][a 3 -6 3 -2ao 2 + a£ 2 ] 2 1 ' 

for the new inverse scattering length and effective range. Here 

£ = K = T^M«M K ( M *\ )\m v {2M k - 2M,)]- 1 / 2 = 17.409M^ 

* =# = -S^-(M^k) 4[(2M - - 2M * r ° /2 - - 6 - 064M - 1 <42) 



at k — 0, and the primes serve to distinguish these scattering lengths and ranges from those appearing in Eq. i|24[l . 
to which they reduce in the no annihilation limit £ — > oo. The numerical values of £ and £' have been calculated 
for meson masses M w = 140 MeV, Mk = 496 MeV, M K * = 892 MeV and coupling constant a s = 2.9. Using these 
together with the values for a and b given in Table [I] one obtains 

a Q = 4.281M^ 1 ; r' a = 3. 2031/^ (isoscalar) (43) 

as predictions for the KK isoscalar scattering length and effective range in the presence of absorption. The corre- 
sponding values for the isovector channel are still 

a\ = -(WOIM^ 1 ; r' x = 3.702M^ (isovector) (44) 

since there is no annihilation contribution from the 1 = 0, 1 = 1 partial wave to the KK scattering amplitude due to 
the identity of the outgoing pions in the isospin representation. Higher partial waves cannot contribute in any event 
due to the contact nature of the annihilation interaction. In fact the ao(980) — ► 2tt decay channel is forbidden by 
G-parity. Thus the isovector scattering length remains real if, as here, the irr] decays are ignored. 
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2. Bound states for I — 

The bound states of Eq. I|33[) are determined |33T] by the complex root(s) of f(k) in Eq. 1351) . with Imfc < 0, i.e. 

(((k)+ik)(k~ia)+i{b 2 -a 2 ) = (45) 
Let k = p — in, k > be such a root. Then the total mass M and decay width T of a bound KK pair at rest reads 

P = M-ir = 2M fe + (p ~ j )2 (46) 

This root can only be found numerically as will be done shortly jH^. However, since the magnitude of the absorptive 
coupling M^\cq\ ~ &wMk/£ ~ 0.7 is quite small relative to the isoscalar potential strength, a perturbative solution 
in inverse powers of C is useful. Replace £ by its value — i£ at k — (more correctly k w — in but the difference is 
quantitatively insignificant). Then Eq. 1451) reduces to a quadratic equation for k. The relevant root is 

b 2 - k 2 

k = — in + 0{C 2 ) ~ (0.1255 - 0.1936i)M K (47) 

where k = ia = —in = — 0.1936-Mr-i is the root leading to the bound state in the absence of absorption, and 
b = 1.491Mk-. Working to the same order in (i.e. neglecting the real part of k), one finds that 

Po « 2M K - XT " 2lK M c ( 48 ) 

Since 47r , !/> 2 (0) = 2k(6 2 — n 2 ) from Eq. (| A. 1 0|> the result in Eq. I48|) is equivalent to treating the annihilation potential 
— CQ<5 3 (r) as a first order perturbation of Po, 

P ^2M K -^---T; r = -^ 2 (0) = — ^ 2 (0) (49) 

Upon recalling that — i^ 1 ~ — J(Po), one sees from this form that the complex solution for Po always lies on the 
second sheet — 2ir < <fi < of the cut complex Po plane for J(Pq). Otherwise the signs of the imaginary parts of 
Po and — J(Po) will differ, and no solution is possible. Taking cognizance of this fact the exact root of Eq. (|45|l and 
resulting value for Po are found numerically to be 

k = (0.1254 - 0.1972i)M K ; P = (981 - 25i)MeV (50) 

using the the previously determined values of a and b shown in Table [I] Notice that the exact value of k lies rather 
close to the perturbation estimate given by Eq. (|47|l 

The exact numerical results for the Bargmann potential that is phase-equivalent to Eq. (|17fl are given in the first 
row of Table [H] We contrast these results with the mass and width calculations based on the Weinstein-Isgur non- 
relativistic potential modelQ, as well as the predictions of Morgan and Pennington|3^|, based on a parametrization 
of the Jost function[5i| with parameters determined from experimental phase shifts. These calculations are also 
compared with the earlier experimental results summarized by the Particle Data Group (PDG)[3^|, as well as the 
recent data from the Fermilab E791 Collaboration 26] . 

In order to extract Bargmann potential parameters that would lead to the mass and width values of either Weinstein- 
Isgur, or Morgan and Pennington, we have turned Eq. I|45|) around and asked what values of a and b produce a root 
k = p — in that gives a prescribed value of the complex total energy in Eq. I|46l) . The answer is 

V + K rxi\ 
a = —K — p (51) 



b = yja? - p(p - + {v + n)(a + n) (52) 

after setting £(p — in) = v — it;. We have treated the experimental data in the same way. These relations allow one 
to establish a direct link between the mass and decay width of the bound KK state, and thence extract the isoscalar 
scattering length and effective range via Eqs. I|4l)|l and 1)41(1 . or Eqs. (|63|l and l|64|) . 

The effective range expansion parameters obtained in this way for the the Weinstein-Isgur calculations, the Morgan- 
Pennington fits, the PDG data, and the Fermilab data, are also shown in Table ITT1 Collectively these results are all 
internally consistent, and suggest that the non-relativistic OBE potential that includes absorption via K* exchange, 
provides a natural description for the properties of the /o(980) scalar meson as a bound KK pair. In particular, 
the prediction of the correct decay width from K* exchange, which is parameter-free and can be treated without 
approximation, adds weight to this interpretation. In this connection it is also interesting to compare with the 
current-algebra prediction[37| of T « 660 MeV for the much larger decay width (600 — 1000 MeV, also dominantly 
into two pions) of the /o(400 — 1200) or a scalar meson, that is considered as a candidate for qq bound state. 
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3. Two-photon decay o//o(980) 

The two-photon decay of the /o(980), considered as a bound kaon-antikaon pair, can be treated by the same 
method. Working in a suitable gauge, the two-photon annihilation potential V 11 is given by the single diagram in 
Fig. |3| with the K K — ► 2ir vertices and pion propagators replaced by the analagous electromagnetic K + K~ — > 27 
("seagull") vertices and photon propagators, after division by two to compensate for crossing. One obtains 

y 77 = -c^5 3 (r) (53) 
Without counterterms, c 77 is formally divergent due to the photon loop. However, the imaginary part is finite, 

T 2 1.27ra\ fri . 



where a — e 2 /4ir « 1/137 is the fine structure constant. The two-photon decay width then follows from the second 
expression in Eq. (|49|l as 



2na 2 a 2 n(b 2 — k 2 ) , , , 

T 77 = 2lmc 2 y K+K ^ = -^-^ K + K _ = L_. > = 5 .59 ke v (55) 



where ipK+K- — ^ < K + K \KK > is the probability amplitude for finding a K + K pair in the /o(980) ground 
state. For a state of good isospin, < K + K~\KK > = l/y/2. One can verify this result independently by noting 
that the coefficient 2-kc? jM"^ = w 77 , where <r 77 gives the low velocity cross section for the annihilation of free, 
charged boson-antiboson pairs of relative velocity v into two photons |38l|. The decay width of the bound state is then 
calculated in the standard way|38j as cr 77 times the flux density, or <t 77 x vip 2 <+K _ , which just reproduces Eq^ (|55|l . 

The numerical value found above for T 77 is an order of magnitude larger than both the estimate given in [39l ] and 
also the currently quoted value of 0.39^013 keV from experiment ^0] . In the molecular picture this discrepancy comes 
about because of the extreme sensitivity of the spatial wave function at the origin to the assumed binding energy of 
the KK pair. In [3^ the author estimates tp(0) « j^M^/ 2 using a gaussian potential fitted to a binding energy of 10 
MeV. This is a factor ~ 3 smaller than our estimate of i/j(0) that corresponds to a binding energy of 18.6 MeV, which 
in turn increases our estimated width by an order of magnitude over that given in|3J| (see Fig. 111(1 . 

On the other hand the 2ir decay width, which is also proportional to ip 2 (0), is correctly reproduced by our wave 
function. Thus it would not seem possible to reproduce both the measured pole position of the /o(980) decaying state 
in the complex plane and at the same time the two photon decay width in the molecular picture without intoducing 
some additional mechanism that reduces the fractional parentage f p = < K + K~\KK > of K + K~ in the /o(980) 
ground state. For, apart from this factor, the branching ratio of electromagnetic to hadronic decays is given by the 
expression 

¥ » ^ - fen^rup - - * uru,. m 

that contains known physical constants if one adopts the SU(3) value g-wKK* = —^/^gp-n-K for the K — > irK* coupling 
constant again. The only less well-determined constant in this expression is g pw7 r- However, the spread between the 
various estimates [53, El of this coupling is not nearly large enough to cause the branching ratio to alter by an order 
of magnitude. This, then leaves f p . Since r 77 /r ~ 10~ 5 experimentally, one would require f p ~ l/vl0 down from 
l/-\/2, suggesting that the assumption of a pure KK molecular ground state of good isospin is oversimplified. 

A competing picture for the structure of /o(980) as a four-quark state has been developed in[5|, which leads to an 
estimate r(/ (980) -> 77) « 0.27 keV H3 that is also (with a particular choice of parameters) not in contradiction with 
experiment. One has therefore to conclude that agreement or otherwise between the calculated and experimental 27 
decay width of /o(980) by itself is too ambiguous for deciding between a molecular versus a q 2 q 2 four-quark structure 
of this meson, and that some mixture of these two extreme pictures is probably indicated. 



4- KK cross sections and complex scattering lengths 



The linear combination of scattering amplitudes of isospin zero and one, f p = (/o + /i)/2 = (S p — 1) /2ik determines 
the s-wave S matrix as S p = (So + S\)/2 for the scattering and reaction cross sections for KK scattering in the 
absence of Coulomb interactions. One has 

a elas = ^\l-S p \ 2 , <j r = ^L(l-\S p \ 2 ) (57) 
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where k is the relative momentum of the colliding pair. The reaction cross section has in turn contributions from the 
(quasi-elastic) "charge transfer" channel, K + K~ *-> K°K°, in addition to the two-pion annihilation channel. One 
finds a r = a ex + a a , with 



— 

k 2 



,Sq — Si ^ 2 



since the annihilation cross section a a is simply the average of the corresponding annihilation cross sections 



rW 



k 2 



(i - \Si\ 2 ) 



(58) 



(59) 



in the separate isospin channels that add incoherently. 

The S-matrix elements of good isospin, Si, appearing in these expressions are given directly by Eq. (|37|l in terms 
of the Jost functions for each isospin channel as 



Si(k) 



(f -k)(k- ia) -{b 2 -a 2 )\(k + ib 



aFj)\k-ib) 



.{£ + k){k + ia) + {b 2 

where the parameters in this expression refer to a specific isospin as indicated. Then 



(60) 



M) - 



4tt 



b 2 + k 2 
4tt£ 



(b 2 -a 2 + ba + k 2 ) 2 + {b-a) 2 ( 2 



a 2 {i + k) 2 + {b 2 
b 2 - 



a 2 + k(£ + k)) 2 



H£ + k) 2 + {b 2 - a 2 + k(£ + k)) 2 



(61) 



Note that ~ 1/As at small k as befits a true absorption cross section. One can confirm the expression for the 
annihilation cross section by considering the flux loss in each isospin channel out of the radial wave function 
ip(k,r)/ f(—k) for scattering, into an infinitesimal sphere located at the origin. Employing the boundary conditions 
<p(k,r) — ► 1, <p'(k,r) — > C as r ~ * used to obtain the Jost function f(k) of Eq. iJSSJ, this flux loss is given without 
approximation by 8tt/ {Mr |/(— ^)| 2 )- This is identical with the right hand side of the second of Eqs. (|61|l after division 
by the incident flux v = 2k/ Mk ■ 

Finally, taking the first option mentioned above Eq. I|rt7l) and expanding expanding k cot <5 i(k), where Si is the 
phase of Si(k) = exp(2i5i), one has 



k cot 6i(k) = -l/a/ f + 



^r e / f k 2 + 



where the complex scattering length and effective range are given by the much simpler expressions 



and 



eff _ 



+ ba — i(b — a)£ 



2(b-a)(£ + ia) 2 



[b 2 - a 2 + ba - i(b - a)£)] 5 



(62) 



(63) 



(64) 



Values of these quantities in the 1 = and 1=1 channels are listed in Table II I II Note that for / = 1 , these 
expressions reduce to those in Eq. 1)25(1 . The Re(aJ ) of course coincide with the previous results, Eq. I)4()[l. for the 
scattering lengths when annihilation is present. From Eqs. (|57|) to (|59ll the zero momentum elastic, charge transfer, 
and annihilation cross sections are given by 



= 47rk 



eff\2 



= 47rk 



e//|2 



4-7T 



-Im{a e J s ) 



(65) 



in the effective range approximation in terms of a^/* 
length in the absence of Coulomb interactions, and a%ff 



eff , 
V + 



eff 
„eff 



)/2, the common K + K or K°K° scattering 
)/2- 



One can show from Eq. (|63|l that the sign of Ke(a e /^) is determined by the sign of the "bare" scattering length 
ai given by Eq. Q25J) in the absence of absorption, and hence the sign of the parameter a of the Bargmann potential 
that controls the existence (a < 0) or not (a > 0) of a bound state in channel /; Im(a/^) is independent of this sign. 
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Fig. 0] illustrates how the real and imaginary parts of a e J ^ would vary should this parameter vary (i.e. the binding 
energy vary) in the isoscalar channel. The remaining parameters for this illustration have been kept the same as in 
Table [I] However, the effects for small variations in the binding energy from real to virtual around zero are not large: 
for example by simply switching the sign of a to +0.9136Mk in Table [I] (thereby also unbinding the isoscalar channel), 
causes o~ e ias to decrease by ~ 20% from 106 to 86 mb while o~ a remains unchanged. 

Calculations of the total elastic, reaction, and v times the annihilation cross sections based on Eqs. (|57|l and Ij58(l 
are shown in Fig.[S] However, as there are as yet no reported measurements of KK scattering with which to compare 
any of these cross sections, we have attempted to use detailed balance to derive the tt + tt~ — + KK cross section, for 
which measurements do exist J2?|. In this case the latter cross section can be exactly related to the 1 = absorption 
cross section for KK — > tttt defined by Eq. H59|) as follows (see Eq. (|A.22|l ) 

pV (7r+7r- - KK) = 2 -k 2 a^ = ^(1 - rg) (66) 

where ryo = |So(A;)| is the KK inelasticity, and where p and k are the relative momenta of the pion or kaon pair of the 
same total CM energy. 

The two curves marked (a) + (b) in Fig. (the upper one includes Coulomb distortion^?] via the factor 
|/i \k,d)\~ 2 ~ 1-KT]\ex'p(2-Kr]) — 1] _1 taken from Eq. (| A. 1 9|) of the Appendix for kd ~ uad <C 1) show the in- 
ferred results for <7 a (7r + 7r~ — > KK), where they are compared with tt + tt~ scattering data[23, |28j, and the fit, curve 
(c), to these data using a parametrized tttt scatterin g am plitude with parameters extracted from these experiments. 
The extraction procedure is documented in detail in[27|]. A coupled channel treatment of the tttt-KK system leads 
to similar results pol l23l l45j (especially |4^) shown by curves (d) and (e). 

Fig. clearly shows that the strong enhancement of the cross section near threshold is not adequately reproduced 
by the pure molecular picture. The same remark holds for the calculated tttt inelasticity, discussed below, that is 
correspondingly somewhat too large near threshold ( Fig.0). 

The pure molecular picture agrees reasonably well with the result obtained by that coupled channel calculation of 
45] which omits the bare scalar meson state /g (curve (d)). This numerical finding justifies the assumptions made 
in order to obtain a simple analytical model. In order to obtain a good agreement with the data, the introduction of 
an explicit bare scalar meson state is necessary (curve (e)), however. 

The parametrized tttt scattering amplitudes as taken from experiment (2?J (curve (c)) assume the existence of a 
pole on the second Riemann sheet just below the KK threshold at [997 ± 6 — (27 ± 8)i)] MeV in the tttt scattering 
amplitude. The experimentally inferred value of this pole lies reasonably close to the value of the pole on the second 
sheet for a decaying /o(980) bound state given in Table [HJ that has been calculated from the very different points of 
view presented here and elsewhere Q- 

A further experimental check on these calculations is provided by comparing with the data [29( for the tttt isoscalar 
inelasticity parameter ^(tttt). Since we have assumed that only coupling between the tttt and KK channels is 
important, this paramater equals the KK inelasticity parameter r)o(KK) = \So(k)\ that may be calculated directly 
from Eq. (|60|l . The resulting prediction for t]q(tttt) is given in Fig. [7| For comparison we also show the behavior 
of 7]o(tttt) that follows from a fully relativistic treatment of the pion loop integral under Pauli-Villars regularization 
with cutoff Apv = 4 GeV (the result is rather insensitive to the actual choice of cutoff). There is no difference in 
the two calculations near threshold. However, the relativistic version becomes more elastic than the non-relativistic 
calculation with increasing CM energy. 

V. KAONIUM 

Besides the isoscalar KK bound system, the K + K~ component of this state can bind to form a kaonium atom 
in the long range Coulomb potential. In the following we construct the Jost function that determines the bound 
states and decay widths of kaonium, that may then be calculated without approximation. The hamiltonian now 
also includes the Coulomb field V c = aQQ/r, (Q,Q are particle, antiparticle charge operators) that breaks isospin. 
However, due to the very disparate length scales ~ MZ and ~ (aMjj)" 1 over which the strong and much weaker 
Coulomb interactions operate, isospin breaking is essentially confined to that region of space r > d outside the range 
of the strong potential. In this region the wave function reads 

^ P \K + K- > +cp n \K K > ), r>d (67) 

where ip p and tp n describe the s-wave radial motion of K + K~ and K°K° separately in the external region, 

d 2 

("T~2 + k 2 - M K V c )ip p = 0, r > d 



13 



d 2 

(j2 + fc2 )^« = °> T > d ( 68 ) 

Inside of this range we can ignore isospin breakin g en tirely due to the dominance of the strong isospin conserving 
interaction (no internal mixing approximation |42l |43j) and work with states of good isospin. To obtain the Jost 
function we therefore need to join the linear combinations — (cp p ± <p n )/ \[2 corresponding to I = and 1 in the 
external region r > d onto the known solutions for the complex potentials in Eq. (|33[) in the internal region r < d. 
The latter are given by tp(k,r) of Eq. (|A.5|I with the appropriate values of a, b and £ taken from Tabled for each 
isospin channel. 

For r > d the solution describing the scattering of a K + K~~ pair is given by the standard expression 

<P P (r) = ftH^r) - S{k)e 2 ^ f^\k,r) (69) 

The fc \k,r) and fc + \k,r) = f* 1 " (fc, r) represent incoming (outgoing) waves in the Coulomb field, S(k) is the 
scattering matrix of non-Coulombic origin, a = argT(l + irf) the Coulomb phase shift and ir\ = [uxfik the Coulomb 

parameter for attractive interactions; fi = Mx/2 is the reduced mass. The /c (fc, r) are proportional to Whittaker 
functions. One has 461 



ft\k,r) = e^/ 2 W imh (2ikr) 

— '2 

— jl - ikr - 2p,ar[ln2ikr + ip(l - irf) + 2-/ - 1] H j, kr -> (70) 



r(i-wj) 

where the small kr behavior of fc (k, r) involves [48) the standard "^-function, and 7 = 0.57721 is the Euler constant. 
The S'-matrix in Eq. I|69l) is determined by the logarithmic derivative C p = (p'/ip p of tp p at r = d. In detail, 

S {k) = e- 2 ^(^M (71) 

where Q c = fc '//i ' is the logarithmic derivative of the incoming Coulomb wave at r = d, and p = arg/i \k, d). 
Thus the Jost function describing kaonium is proportional to £ p — £ c . To find it [43} we determine the value of £ p by 
matching tp p ± ip n onto the internal solutions ipo and ipi of good isospin, where (p n ~ exp(ikr) is a pure outgoing 
wave [54). Due the large difference in space behavior of the strong (br ^> 1) and Coulombic (kr -C 1) wave functions 
near r ~ d for momenta fc ~ [ia of interest for kaonium, the <p(k, r) have already reached their asymptotic behavior, 

while the fc (k, r) are still given by their small kr approximation. The matching of value and derivative at r = d is 
therefore particularly simple to perform. One finds 

1 »" ^ (72, 



Cp 2 V Co £ x , 

where the logarithmic derivatives (j = tp'j/(pi of the internal wave functions carrying good isospin / are 

'/(fc)e 4M + /(-fc)e- lM 



" V/(fc)e* M -/(-fc)e- 4fcd // 1 J 

The eigenvalue condition £ p = Q that determines complex energies of kaonium then reads 

1 - ~Ce(i + h = (74) 
2 Co ft 

This equation contains no free parameters. The Jost functions /(fc) in the C/ all have pre-determined constants for 
each isospin channel as given in Table |H and come from Eq. I|35|l if annihilation is included, or Eq. 1221) if it is not. 

It is useful to express the roots of Eq. I|74|l in the form k\ = — i\/j,a since then A -1 — > 1, 2, 3, ... for a pure Coulomb 
field, and also to recognise that the logarithmic derivatives ft can be replaced by their zero momentum values with 
impunity because fc is always O(fia) for kaonium. Then 

^«-(af f -d)+0(l/C 2 ) (75) 
ft 

in terms of the complex scattering length of Eq. (|63|l . 
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Using these forms for the (j , the eigenvalue condition i|74[l can be reformulated in terms of the complex scattering 
length ay' for K + K~ scattering as 



l(tf' + a?') = (l-^)d (76) 
This leads to an eigenvalue equation of the Kudryavtsev-Popov (KP) tvpepn). 

l/a e p Sf = 2fia ip(l - 1/A) + ^A + ln(2A/zacZ) + 2 7 (77) 

The logarithmic derivative Cc(k) of the Coulomb wave function has been calculated at k\ — —iXfia from Eq. I|70jl . 
The resulting complex energy spectrum for kaonium is then given by 

E x - lr A = -\\ 2 M K c? (78) 

The solutions of the KP equation are dominated by the -0-function on the right hand side that has poles at A = 1/n, 
for integer n. These poles just reproduce the pure Coulomb spectrum of kaonium. The strong interaction changes 
this behavior by moving A off these poles to A = (n + <5) _1 , where 5 is complex and not necessarily small. In Table llVl 
we list a series of exact numerical eigenvalues A of the KP equation, and the associated complex energy spectrum for 
kaonium. These computations use the value 



a P 



l ff = (1.940 - 1.199i)M^ x = (0.77 - 0.48z)fm (79) 



for the complex scattering length calculated from Eq. I|63|) plus the data given in Table IIIII Notice that the ao (980) 
admixture introduced via the isovector scattering length a\ only contributes to the real part of the K + K~ effective 
scattering length that enters into the calculation of the level shifts and decay widths of kaonium. The imaginary 
part is entirely determined by the isoscalar piece of the strong interactions. However, due to the non-linear nature 
of the eigenvalue problem, the isovector contribution nevertheless influences both the real and imaginary parts of the 
resulting energy spectrum. Only in the perturbative limit (to be described below) does the ao(980) mode contribute 
solely to the level shifts. 

For comparison we have also used the isospin zero and one scattering lengths a,Q = (0.86 — 1.82i) fm and a\ = 
(0.06 — 0.77i) fm, based on a recent investigation into the effects of a high lying scalar isoscalar quark antiquark 
configuration (which we refer to as /g) mixed into the KK molecular ground state. Here the effect of the 7r7y-channel 
is included too. The K + K~ scattering length then becomes 

a e Jf = (1.156 - 3.268i)M^ x = (0.46 - 1.30z)fm (80) 

instead of Eq. (J75J . 

The predictions for kaonium are given in Tables HH and E| while Table EH summarizes the decrease in lifetime 
of the ground state as the complexity of the K + K~ interaction is increased. In particular, Table VI shows that 
the lifetime of kaonium is reduced further by about a factor of three once annihilation in the isovector channel is 
introduced via the irr/ decay of the eto(980). These results are qualitatively similar to those given by other approaches 
to the problem of kaonium formation and decay0|- One notes from Tables llVl and E] that ReA < 1/n, corresponding 
to a positive shift of the Coulomb level when the strong interaction supports a bound state (the sign of the real part 
of the scattering length determines the sign of the level shift). In addition to causing repulsive shifts of the Coulomb 
levels the strong field introduces a decay width due to K + K~ — >■ 2ir annihilation with a lifetime of ~ 10~ 18 sec (the 
two-photon annihilation lifetime at 0(M^ a -5 ) ~ 10 -13 sec, or one in 10 5 decays, is too long to compete). It also 
introduces an extra node into all the eigenfunctions of kaonium. This means that the kaonium levels represent excited 
states in the combined Coulomb plus strong fields of the KK system. As will be discussed in more detail below, 
the occurrence of this node is responsible for making the kaonium lifetime two orders of magnitude shorter than the 
~ 10 -16 sec it would have been had the K + K~ pair only been bound by a Coulomb field 55] . 

Let us briefly indicate how these features come about. The appearance of the additional radial node is required by 
the fact that the Is state of relative motion has been usurped by the KK ground state in the strong potential. So the 
ground and excited levels of kaonium have to have a 2s, 3s... character with one additional node. This node is always 
positioned in the close proximity of the K + K~ scattering length, independent of the excitation energy. It is real or 
virtual depending on whether the strong interaction supports a real (a p > 0) or a virtual (a p < 0) bound state. This 
comes about as follows. The kaonium levels are all weakly bound at k ~ [ia so that their radial wave functions in 
the inner region all become nearly identical with the zero energy wave functions of KK in the strong potential alone. 
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These functions [see Eq. (|A.15|) ] all have a node at r = a p in the effective range approximation. The remaining nodes 
are supplied by the Whittaker function 1 (2\fiar) in the external region, r > d. They are almost identical with 

those of a pure Coulomb field bound state wave function where A = 1/n. 

The eigenfunction for the lowest state of kaonium is illustrated in Fig. [5] (with absorption suppressed). It is clear 
from this figure that the actual choice of the cutoff d is immaterial if taken anywhere in the region where the wave 
function behaves linearly. A closed form for the bound state eigenfunctions (fii/\ of kaonium labelled by the quantum 
number A -1 is given in Eq. (|A.13|I in the Appendix. From that one calculates that the gradient at the origin is 
enhanced by the factor b/a in both isospin channels over that of the pure Coulomb state. Physically this comes 
about because the strong potentials draw the K + K~ pair much closer together, thereby increasing the probability 
density of finding them at the origin by (b 2 /a 2 )i (notice b > \a\ from Eq. (|26|l ) in each isospin component, where they 
annihilate. In the isoscalar channel, which is strongly attractive, the effect is especially marked, (b 2 /a 2 )o ~ 60. On 
the other hand the isovector channel only supports a virtual state. Hence the concentration of probability density at 
the origin in this channel is much weaker, {b 2 / a 2 )\ ~ 3. 

The effect of this factor in the isoscalar channel is quite dramatic. It increases the decay width of kaonium by two 
orders of magnitude (with a corresponding decrease in lifetime) over what it would have been in the absence of the 
strong potentials Vj. One can illustrate this qualitatively by appealing to the perturbative solution [30ll3l| of the KP 
equation for the complex energy shift. For the ground state this reads 

AE l -T « = \u 2 K o?a v - \*fo?C X [|(£)o] (81) 

that equals either (0.187 — 0.116i) keV, or (0.187 — 0.164?) keV depending on whether the value (79) or its low 
absorption version £ 3> 1 taken from Eq. I|63|) , 

a eff „ K + Ot) _ *_ { #_ )o = (L940 _ L703)M -1 (82) 

is used for the effective scattering length. The aj are the bare scattering lengths of Eq. (|25[) as before. The last 
approximation places the role of the isoscalar probability density admixture (b 2 /a 2 )o in the kaonium ground state 
explicity in evidence. 

Neither approximation is very accurate in comparison with the exact solution of (0.181 — 0.107z) keV, but can be 
used to illustrate the effect of strong interactions decreasing the lifetime of kaonium against what it would have been 
in a pure Coulomb field. One sees this as follows. The decay width in Eq. (|81|) is directly related to the "annihilation 
at rest" cross section 

47T ff 47T 1 r 6 2 £ , 47T r l / & 2 , 

a a « - T Im(a/f) = - - - ^ - ^ « ^[5(^)0] (83) 

as r = v<fi 2 (0)a a , upon multiplication by the flux factor v x 2 (O) where 2 (O) = M^ct 3 /8ir gives the probability per 
unit volume for finding the K + K~ pair at the origin when bound by the Coulomb field alone; v = k/n is their relative 
velocity. On the other hand when the strong potential is absent, the scattering lengths a^^ — ► — i/£ and — > are 
pure imaginary and zero respectively, so that a e J^ — > — i/2£ — — 0.028iM^ , and the corresponding free annihilation 
cross section becomes 

47T 1 

Vajree = ^t(-) ~ (T a /9, ~ 42 to 60 (84) 

if we use the unexpanded or expanded version of a a in Eq. (J53J). The factor reflects the modification of the kaonium 
ground state due to the strong interactions. 

In this case there is no energy shift at all, and the decay half-width is ^Tf ree = 0.0028 keV, leading to a lifetime 
of Tf ree — 1.2 x 10~ 16 sec. So the K + K~ — > 2ir strong potential assisted annihilation lifetime drops by roughly the 
factor 8 from 1.2 x 10~ 16 sec that it would have been in a pure Coulomb field, to 3.2 x 10~ 18 sec as given by the 
exact solution of the KP equation. Table IVll summarizes the change in lifetime of the kaonium ground state as the 
complexity of the K + K~ interaction changes from assuming free annihilation at rest in a Coulomb field, to including 
one-boson exchange, to adding /q mixing, for determining the scattering length. 

We note in passing from the center expression for a a in Eq. I|83|l that if the bound state occurs at zero energy 
(a = 0), the corresponding enhancement of the annihilation width is maximal. Also, while generally true of the exact 
solutions of the KP equation, the perturbative solution given in Eq. Ij81|l makes it evident that the sign of the level 
shift goes along with the sign of the scattering length, which is positive for bound and negative for unbound or virtual 
states. 
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We close this section with a comment on to what extent the results for kaonium might be considered as " parameter- 
free" . Apart from specifying a Lagrangian with the attendant assumption of SU (3) symmetry that fixes the various 
coupling constants in terms of g p7r7r , that we in turn fixed from the KSRF relation, the only arbitrary parameter in 
these calculations was the form factor cutoff A introduced in Section IIIA. The choice A = 4 GeV, that was also used 
in Ref. |2f|. led to calculated scattering lengths and effective ranges for the KK system, that were then duplicated 
by constructing an equivalent Bargmann potential that depended on two parameters a and b. This in turn gave 
predictions for the mass and width of a bound KK molecular state in reasonable accord with the known experimental 
mass and width of the /o(980) scalar meson. However, as already demonstrated in Table ITT1 one can turn the argument 
around without materially affecting any of our conclusions for kaonium, by asserting that the Bargmann parameters 
a and b are actually pre-determined from experiment in the strong interaction sector, thus leaving the kaonium sector 
without any free parameters at all. 

Similar remarks hold true for the calculated cross sections shown in Figs. |S] and |SJ Once a and b have been fixed 
from experiment or otherwise, no further adjustable parameters enter into these calculations. The actual values of a 
and b as determined by any of the five alternatives shown in Table ITT1 do not in fact differ very much. 

VI. DISCUSSION 

In this investigation we have shown that one meson exchange potentials derived from a standard SU(3)v x SU(3)a 
Lagrangian in the non-relativistic limit are sufficient to bind KK into a kaonic molecule with a mass and decay 
width that closely match the experimental values of the /o(980) isoscalar meson. We have also shown that the same 
potentials, in combination with detailed balance, lead to a prediction of the ir + n~ — > KK cross section that is too 
low especially near threshold. 

We have further demonstrated that there are specific, experimentally observable effects in the energy level shifts 
and decay properties of the kaonium spectrum, in the presence of a KK molecular ground state. The existence of 
such a bound KK pair would be reflected in positive energy level shifts in the kaonium Coulomb spectrum (due to 
level repulsion), plus enhanced K + K~ — > 2ir decay widths. On the other hand, unbound or virtually bound (and 
therefore weakly correlated) KK pairs produce negative energy level shifts and a weak enhancement of the decay 
widths. Since the sign of the level shifts and the enhancement of the decay widths go together with the presence 
or absence of a strongly correlated KK pair in the ground state, measurements of the level shifts and decay widths 
in kaonium could provide another probe for investigating to what extent the Jo (980) meson ground state has a KK 
molecular component. From the experimental point of view it is important to stress that, as for protonium|3lj which 
has a very similar decay width, the Is — 2s level spacing of ~ 5.1keV in kaonium is many times the estimated decay 
widths < 200 eV, so that these levels remain distinguishable in the presence of the 27r decay channel. 

Notwithstanding the results reported in (3^, we again stress the conflict that exists between the pure molecular 
description of the /o(980) ground state and the calculated 27 decay width. In 39] a binding energy of 10 MeV for 
the kaonium molecule was assumed, that leads to a width of ~ 0.6 keV. However, most recent data show a resonance 
energy that corresponds to a pole in the complex plane, which translates into a binding energy of 18 to 20 MeV for 
the molecular potential. This increase in binding energy enhances the value of the spatial wave function at the origin 
by a factor of three, see Fig. 1111 One way of understanding the origin of this discrepancy is to generalise the molecular 
picture of the /o(980) to a superposition of a KK component and a hard component. In the Jiilich model [20( a f' Q 
state has been considered in addition to the pure molecule in order to reproduce the large inelasticities of the 7T7t phase 
shifts in the vicinity of the KK threshold. The /q state may be interpreted as either a qq or as a q 2 q 2 configuration. 
As pointed out by Achasovj^ and confirmed by current algebrajS^, the large decay width of the former into two pions 
may disqualify a qq interpretation of the additional component f . The effect of the f' a component in the /o(980) 
meson introduced in the Jiilich model is to reduce the lifetime of kaonium from 3.2 x 10 -18 sec to 1.1 x 10 -18 sec. 
The presently available data concerning the inelasticities of the ttt: phase shifts, and in particular the 7r + 7r~ — > KK 
cross sections, do not allow one to disentangle the components of the /o(980) structure in a model independent way. 
High resolution data for two kaon production near to the threshold will be undertaken at COSY in the near future. 
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APPENDIX 
1. Isospin phase conventions 

In order to construct two-meson states of good isospin using conventional Clebsch-Gordan coefficients the following 
phase conventions are useful. 

If one chooses |3^| as basis states 1 7r ± >= =p|7r, ±1 >, and \ir° >= \ir, > for charged and neutral pions, the two-pion 
state of good isospin \I = 0, I3 = > is, for example, 

IHo >= -y|[k+(l)7r-(2) > +|7r-(l)7r+(2) > +|7r°(l)7r°(2) > ] (A.l) 

Since the pions are identical bosons in the isospin basis, one has to include a space part of the two-pion wave function 
that is symmetric in particle label interchange to accompany \I = 0, 13 = >. The normalization of the symmetrized 
product of the two-pion space wave function introduces an additional factor n s = 1/^/2. 

For the kaon and anti-kaon isospin doublet and anti-doublet the analogous convention is, in matrix form, 

Consequently the isoscalar and isovector combinations \I = 0, I3 = > and \I = 1, 13 = Q > are 

\(KK) >= -±=[\K+(l)K-{2) > + \K°(1)K°(2) >} 

\{KK) X >= -±=[\K+{l)K-{2) > -\K°(1)K°(2) >] (A.2) 
These relations lead immediately to the values of the isospin coefficients Cj and n s quoted below Eq. I|13[) in the text. 

2. Wave functions for the Bargmann potentials 

We briefly summarize some results of Bargmann for constructing phase equivalent potentials. As shown in 24] the 
assumed form Eq. (1221 for the Jost function gives rise to a family of potentials that lead to this form for f(k): the 
simplest one of this family is given by Eq. (|23|) . Let ip(r) — r -1 ^^), and replace the OBE potentials by this form for 
Vi. Then 

tp" + (fc 2 - M K Vi)ip = (A.3) 

for K K s-wave scattering. The two linearly independent solutions /(±/c, r) for this equation that are irregular at the 
origin arc known in closed form by construction |24l l33| 

_ ikr k — ib[bsinh(br) + acosh(6r)][6cosh(6r) + a sinh(6r)] _1 

J(k,r) =e — 

k — 10 

^e- tkr , r->oo (A.4) 
A general solution of Eq. (|A.3J) can be written as a linear combination of the /(±/c, r), 

V(k, r) = JL[/(-fc)/(fc, r ) - f(k)f(-k, r)] (A.5) 

where the coefficient /(fc) is the Jost function. It can be shown quite generally |33j that /(fc) is determined by the 
Wronskian of f(k,r) and ip(k,r): 

/(fc) = W[f(k, r), <p(k, r)) = \}<f? - /V] (A.6) 

Since W is independent of r this can be evaluated at any convenient point, in particular at the origin. The value of 
/(fc) is determined by how ip behaves at the origin. If (p ~ r, then /(fc) = /(fc,0) = (k — ia)/(k — ib) as in Eq. 
This is the usual case. On the other hand if ip does not vanish at the origin, but instead possess a finite logarithmic 
derivative there, then /(fc) is given by Eq. 1)35(1. 
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In the first case one obtains the usual scattering solution of Eq. (|A.3|) that vanishes at the origin as 

sinfcr b 2 — a 2 1 fcsinh(br) cos(fcr) — bcosh(br) sin(fcr) 
^ ,T ' = ~~k~ + b 2 + k 2 k asmh{br)+bcosh{br) ( ' ' 

If a = — k is negative, then f(k) has a zero in the lower half of the complex fc-plane, corresponding to a bound state 
with binding energy e — k 2 /Mk- The associated bound state wave function is found by setting a — —k and k = —in 
in Eq. (fO)l to find 

ocosh(or) — /tsinh(or) 

The normalization of this bound state is also available from a knowledge of the Jost function: 

Vc-fa.r) - - = (A.9) 

since /(fc) = [k + in)/{k — ib). The normalized bound state is then given by 

^(r) = (V^ryWMb 2 - k 2 ) u ^ Hbr) . Uh , e~ Kr (A.10) 

Ocosh(Or) — /csinh(or) 

In Fig. 1111 we compare the radial eigenfunction as obtained by numerically solving for the KK isoscalar bound 
state in the one boson exchange potential of Eq. 1)170. with the analytical solution given by Eq. I|A.10|I for the phase 
equivalent Bargmann potential. Except in the vicinity of the origin, where the numerical solution is underestimated 
by a factor ~ 1.2, the two wave functions are essentially identical. Also shown is the wave function at 10 MeV binding 
in a gaussian potential that was used in |39| to estimate ip(0). 

Finally in Fig. Owe illustrate the variation of the spatial wave function ip(0) of Eq. IjA.lOf) as a function of its 
binding energy. This plot was obtained by using the calculated variation of eto, r$ and e with cutoff as shown in Fig. [2 
for the one boson exchange potential to obtain the parameters a — — k and b at each binding energy. 

The inclusion of the annihilation potential changes the solution for the radial eigenfunction from Lp(—in, r) to the 
complex solution 

(p(k r ,r) = f(k r ,r)/f(k r ,0) (A.ll) 

where k r = (0.1254 — 0.1972i)Mk is the root of the Jost function in the presence of annihilation as given by Eq. JSUJ. 
One readily verifies that <p(k r ,0) — 1 and (p'(k r ,Q) = C = ~ as required. A comparative plot of the normalized 
versions of the moduli of (p(k r ,r) and ip(—iK,r) is given in Fig. El Since the derivative of <p(k r ,r) at the origin is 
purely imaginary, its modulus has a zero derivative there. This is clearly visible in the inset in Fig.|3 that also shows 
how \if(k ri r)\ "heals" within 0.1 fm away from the origin from the effects of the annihilation potential to rejoin the 
real solution ip(—iK,r). Either eigenfunction leads to a rms molecular radius of \/ < r 2 > « 1.62 fm, (i.e. about two 
proton charge radii). 

The momentum content of the normalized spatial wave functions 4>{r) = {\fA^r)^ 1 ip{~in 1 r) and ip(r) = 
{^/~Aftr)~ ip(k r ,r) is given by their fourier transforms i/j(p) and tp(p). These are shown in Fig. 1101 for relative me- 
son momenta up to 1 GeV. As expected from Fig. the effect of the annihilation potential is to redistribute the 
momentum to higher values, while at the same time leaving the characteristic value of < p 2 > at ~ 0.08 GeV 2 for the 
relative meson momentum squared in the molecule. 



3. Matching of isospin amplitudes 

The wave functions of kaonium in the inner, or isospin conserving region are given by Eq. (|A. 5|) . These solutions 
either vanish at the origin, or possess a given logarithmic there, depending on the form of f(k) that is chosen. Matching 
either form onto the linear combinations ip ± ip n of the outside solutions at r — d, one finds the consistency condition 

( Co - Cp Co - Cn \( \ _ n 

Ul-Cp -Cl+Gj \<Pn) 

where £ n = tp' n /(p n . The other logarithmic derivatives have been defined in the main text. At momenta k ~ fia of 
relevance for kaonium, £„ = ik is much smaller than any of these so we simply drop it and find 

^~\{^ + ^); ¥>«M>«-1 + Cp/Co = 1-Cp/Ci (A.12) 

Sp A 1,0 SI 
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The eigenvalue condition £ p = £ c determines the values of k = —in, which are pure imaginary in the absence of 
absorption, for the bound states of kaonium, and at the same time removes the outgoing wave piece of ip p from 
Eq. (|69|1 . Then the matching conditions show that the kaonium eigenfunctions that join smoothly onto the pure 
decaying state /i \~inr) ~ i(2nr) outside the strong potential are given by a linear combination of isospin 

zero and one internal functions as 



<Pi/\(r) = TjCcfc 



1 tp (-in,r) 1 tpi(-in,r) 



Co ¥>o(—iK, d) Ci ^(-w. d) 



< d 



= f c (-iK,r), r>d (A.13) 

where Cc = Ce( — an d fc = fc (—iK,d). We have labelled these eigenfunctions with the "quantum number" 
A -1 = fia/n that reverts to an integer for a pure Coulomb field. One verifies immediately that the matching of value 
at r = d is guaranteed by the eigenvalue condition (|74|) or (|A.12I) . while the derivatives match identically from the 
definition of the logarithmic derivatives. 

These eigenfunctions develop a common additional node at the K + K~ scattering length, independent of the energy 
of the kaonium level. Consider the case where the scattering length is real. Then this happens because k ~ /J,a <C Mk 
is very small so that the behavior of the <fi{— in, r) differ but little from their zero energy behavior. From Eq. (jA. T|> 
this is 

, . b — a t (b + a) sinh(6r) ^ ba b sinh(for) + a cosh(6r) | 
b 2 \b cosh(6r) + a sinh(6r) b — ab cosh(6r) + a sinh(6r) / 
~ r, br <C 1 

~-(r-aj), 6r>l (A.14) 
where a and b refer to channel I. Inserting this information into Eq. (|A. 13|> . one finds that 

¥>i/\ (Ccfflp)/c(l - r/dp) (A.15) 

outside the range of the strong interaction. 

We return to the neglect of the K°-K ± mass difference A in calculating the logarithmic derivative £ n for the K°K° 
channel. Including ( n = —t/2MkA leads to a revised version of the eigenvalue equation (|74^l that now reads 

l-kc + C„)(-^ + -^) + ^ = (A.16) 
z so y 4o4i 

The lowest eigenvalue of this equation is A = 0.9860 + 0.0061? instead of 0.9863 + 0.0082i, an inessential change. 
The reason is the dominating influence of the isoscalar channel for which £ d » 18.0 — 10.5z, that overshadows the 
remaining logarithmic derivatives, which are all 0(1), and suppresses the last term in Eq. I|A.16() . 

4. Coulomb distortion 

It is also useful to document the modification due to the Coulomb interaction on the low-momentum behavior of 
the absorption cross section as given by the second of Eqs. JfJSJ). Using Eqs. (|71|l . (|72l) and J7SJ), one finds that 

a a « Im( a ;fl)Im{Q) » - — , , p ' (A.17) 
fc2 P k \ft\k,d)\ 2 

since 

Jm(Cc) = -k\ft\k,d)\- 2 is the imaginary part of the logarithmic derivative Cc of the incoming Coulomb wave. 
The latter relation is easily found by writing /j (fc, r) = exp(ztr)(G(fc, r) — iF(k,r)) in terms of the regular and 
irregular Coulomb solutions F(k, r) ~ sin(fcr — r)ln.2kr + a), G(k, r) ~ cos(kr — r]ln2kr + a), so that 

(GG + FF' -ik\ ,. , a . 

C ° = ( G 2 + F 2 L, (A ' 18) 

after using the Wronskian relation (GF' — G'F) — k that holds for all r. For kd ~ [iad << 1 one finds from Eq. I|70|) 
that 



l/3-'ft ">l 2 = (^) [(> - + 2,-1))% ( ? |^) 



(e 2 ^ - 1). 

(-2^)' 5 = -W* (A.19) 
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in agreement with Rcf. 47 



Detailed balance 



If one ignores the real one-boson exchange potential Vj in Eq. 1|33[1 . then the Jost function for KK scattering due 
to the annihilation part alone is simply f(k) = £ + ik. This gives rise to a s-wave 1 = scattering length of 



Ana (KK) 



-M K c 2 



(A.20) 



after using the definition (I34|) for (. However, since irir annihilation into KK is described by the same diagram in 
Fig. |3 with the pion and kaon lines interchanged, the 1 = 0, tttt scattering length due to this process alone is also given 
to the same approximation by the expression as above, but with M n and Mr interchanged in the product MkCq, i.e. 



47ra (7T7r) = -[M K ($] Mw ^ Mk 



-M K cl = 4ira {KK) 



(A.21) 



We have used Eq. (|31l) to perform the mass interchange. Here p and k are the relative momenta of a pion or kaon pair 
with the same total CM energy Pq. There is no contribution to the a2(7T7r) scattering length from Fig. [31 in the 1 = 2 
channel, since the intermediate kaon loop can at most carry 1=1. Consequently detailed balance holds between the 
1 = absorption cross section cri°^(7r7r — > KK) and its inverse era (KK — * ^ir) when viewed in the isospin basis 
under the above approximations. 

However the physical tt + tt~ channel is described by the scattering length a(7r + 7r~) = |ao(7T7r) + ^0,2(1^). But 
7m(a2(7T7r)) = 0, so 



£r a (7r + 7T 



4tt (2 . , 
KK) = Im -ao(7T7r) 

P V 3 



provided that Pq > 2Mk- 



(A.22) 



6. Dimensional regularization of J(Po) 
Replace the three-dimensional integral in Eq. (|29|l by its analog in d dimensions [T^. 



1 



or, as d 



M(2M - Po) 



J(Po) = --r M{2M~P Q 

47T 



3 \ l 2 /M + 2M - P 

d/2-1 m 



(4 7r )d/: 



(4^/2 r(d/2) J a 
r[l - d/2] 



dL" 



L 2 /M + 2M- Po 



1/2 



M 
Air 



M\(P -2M)\ 



1/2 



i(0+7T)/2 



< 6 < 2tt JP 



(A.23) 



(A.24) 



on the first sheet of the complex Po plane. Note that J(Po) is pure imaginary and positive on the upper lip of the 
cut (j> = extending from 2M to 00. The analytic continuation to the second sheet through the cut starting at the 
branch point Pq = 2M is obtained by taking — 2tt < <j> < 0. 
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[51] This can also be interpreted as <fi(k, r) behaving discontinuously at the origin, tp(k, 0) = 0, but lim p(k, r) ^ as r — > 
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FIG. 1: One boson exchange potentials using a cutoff of A = 4 GeV. The coupling constant is a s — 2.9 and the masses of 
the exchanged mesons are (M p , M^,M^) = (769, 783, 1019) MeV. The analytical Bargmann potentials that reproduce the same 
scattering lengths and effective ranges as the original one boson exchange potentials are shown as companion curves. 



under the influence of the delta function potential. 
[52] For the values of a and b and the logarithmic derivatives given in Table [I] it turns out that there is only one such root in 

the isoscalar channel, and none in the isovector channel. 
[53] Note: these authors use the opposite convention for the definition of the Jost function to that employed hi p^ . Thus their 

sheet II for complex k corresponds to Im k < in our notation. 
[54] This channel is actually closed and decays exponentially like ~ exp(— y/2M K o Ar) when isospin breaking due to the K°-K ± 

mass difference A ~ 4MeV is not ignored. This is an inessential approximation in the presence of the dominating isoscalar 

channel. Keeping A changes the eigenvalues of Eq. I|770 by at most ~ 10%, see below Eq. IIA.160 , 
[55] The predicted lifetime of kaonium is three orders of magnitude shorter than that for the pionium — > 2ty° decav llal . Apart 

from specific details of the annihilation interaction for pionium, its decay rate is slowed down by the relatively small mass 

difference M w + — M n o ~ 5 MeV that drives the decay as opposed to Mk — M n ~ 356 MeV for kaonium. Also, the mr 

system has no bound state analogous to KK. 



23 




3.0 3.5 4.0 4.5 5.0 

Cutoff A (GeV) 

FIG. 2: Dependence of the K K binding energy e (left hand scale), scattering length ao and effective range ro (right hand 
scale) on the choice of cutoff A used in the form factor for the isosalar one boson exchange potential. 



TABLE I: Parameters of the Bargmann potentials, Eq. 11 2 .'-> I . that give rise to identical scattering lengths and effective ranges 
as the one boson exchange potentials in Eqs. 1171 and 118H . The last column gives the value of the logarithmic derivative that 
determines the strength the annihilation potential. The units are Mk- 



Isospin 


a 


b 


c 


1 = 


-0.1936 


1.491 


-17.409i 


I = 1 


+0.9219 


1.462 


-26.114i 
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FIG. 3: Replacement of the irreducible KK scattering diagram involving the exchange of two K*'s by an equivalent pion loop 
diagram with point vertices. 
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FIG. 4: Real and imaginary parts of the complex scattering length a^f as a function of the parameter a of the isoscalar 
Bargmann potential. The remaining parameters are as in Table The /o(980) ground state becomes unbound as Re(a^^) 
passes through zero from positive to negative values. 



25 



TABLE II: Summary of calculated isoscalar KK scattering lengths a' and effective ranges r' (in units M K ), based on the 
a and b Bargmann potential parameters (in units Mk) that reproduce the complex bound state at Po shown in column two. 
The physical constants that enter into these calculations are: (M n , Mk,Mk* ) = (140, 496, 892) MeV, and a B = gp 7r7r /47r « 2.9, 
using the KSRF value of g pv -n ~ 6. The last four rows list the values of a and b, and the resulting isoscalar scattering lengths 
and effective ranges that would be required to reproduce the theoretical results of Weinstein and Isgur ( WI) , the fits of Morgan 
and Pennington (MP), and the experimental data. 

Source P = M - iT/2 (MeV) a b a' r' () 

This calculation 981 - 25i -0.194 1.491 4.28 3.20 

WIf6] 972 - 16i -0.213 1.170 5.01 2.12 

MPf36] 970 - 24i -0.231 1.367 4.30 2.16 

PDG[35J (974 ±3)- (24±5)z -0.22 ±0.02 1.41 ±0.16 4.35 ± 0.49 2.46 ± 0.35 

E791 [26] (975 ± 3) - (22 ± 2)i -0.21 ±0.02 1.37 ±0.05 4.52 ± 0.20 2.51 ±0.33 



TABLE III: Complex scattering lengths and effective ranges for isospin zero and one scattering channels as calculated from 
Eq. 1631 . as well as from coupled channel calculations with the f' and the 7rr;-channel0| included (in units of M^ 1 ). 





Isospin 


a/ 


iff 


One boson exchange 


1 = 


4.281 - 2.398i 


1.169 -0.178z 




I = 1 


-0.401 


3.702 


±/o mixing ±7T7;-channel 


1 = 


2.165-4.592i 


0.151-0.404i 








I = 1 


0. 147-1. 944i 


-1.463-0. 202i 



TABLE IV: Complex level shifts and lifetimes due to the strong interactions away from the point Coulomb value — ^a 2 Mk /n 2 = 
— 6.607/n 2 keV for the ground and excited s-states of kaonium. Numerical solutions of the Kudryavtsev-Popov equation have 
been generated for the K + K~ scattering length a^^ = (1.940 — 1.199i)M^ 1 given by Eq. 1)79^ . The matching radius is 
d ~ 3M^ . The real part of A lies close the sequence 1/1, 1/2, 1/3, 1/4, • • • for a pure Coulomb spectrum. 

Level A AP-ir/2 (keV) Lifetime (x!0~ 18 sec) 

3 rd 0.2491 ± 0.0005i 0.003 - 0.002i 199 

2 nd 0.3318 ± 0.0009i 0.007 - 0.004i 84 

1 st 0.4965 ± 0.0020i 0.023 - 0.013z 25 

Ground 0.9863 + 0.0079i 0.180 - 0.103z 3.2 



TABLE V: A repeat of the calculations described in Table HVI but for the scattering length d e / f = (1.156 - 3.268i)Mj^ of 
Eq. I|80|l which has been taken from the model of Ref. 45] that includes /q mixing and the ±7rr;-channel. This choice leads to 
lifetimes of about a factor three shorter than those given in Table ITV1 

Level A AE - iV/2 (keV) Lifetime (xl0~ 18 sec) 

3 rd 0.2494 ± 0.0014i 0.002 - 0.005i 71 

2 nd 0.3322 ± 0.0025i 0.005 - O.Olli 30 

1 st 0.4975 ± 0.0056i 0.017 - 0.037i 8.9 

Ground 0.9899 ± 0.0223i 0.137 - 0.291i 1.1 



TABLE VI: A comparison of predicted lifetimes of the kaonium ground state for increasing complexity of the K + K interaction. 





Scattering 


; length (frn) 


t(K + K — » tttv) (sec) 


Free annihilation at rest 


0.00 


- O.Oli 


1.2 x 10" 16 


±One boson exchange 


0.77 


- 0.48i 


3.2 x 10 -18 


+ /o mixing ±7rry-channclf45] 


0.46 


- 1.30« 


1.1 x 10~ 18 




FIG. 5: Calculated total s-wave elastic, reaction and tix the annihilation cross sections as a function of the invariant mass Po 
of the colliding kaon-antikaon pair. The vertical line indicates the threshold value at 2Mk- 
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FIG. 6: The -k + tt~ — > KK cross section given by detailed balance, Eq. Ij66^ . as a function of the TfA" invariant mass Po- 
The two curves designated (a) + (b) give the results with (upper curve) and without (lower curve) Coulomb distortion. Also 
shown are the tv + tv~ — > KK experimental data |27|l (closed circles), and |2S| (crosses). Curve (c)) is a fit to these data using 
the parametrized scattering matrix of |2^| . Curves (d) and (e) correspond to coupled channel calculations of the -k + tt~ — > KK 
cross section with (curve (e)) and without (curve (d)) /q mixing. 
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FIG. 7: The -kit inelasticity parameter Tjoimv) as a function of the CM energy. The experimental data is taken from |2£ 
upper curve shows the behavior of the inelasticity given by a fully relativistic calculation of the pion loop integral. 
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FIG. 8: Kaonium ground state radial wave functions: (a) with and (b) without the meson exchange fields present. In (a) the 
strong field solution, Eq. ljA,130 . is joined smoothly at r — d onto the Whittaker function W 1 ^ x j_(2A/iax) for A = 0.981. The 

node in this solution falls at the K + K~ scattering length a p . In (b) the solution W 1 i(2fiar) = 2/xar exp(— fiar) for a pure 
Coulomb field, A = 1, is plotted for the same range of r. To set the scale, note that the Bohr radius of kaonium would fall at 
~ 300 units to the right of this figure. 
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FIG. 9: A plot of (the moduli of) the normalized radial wave functions (p(k r ,r) and ip{— in, r) that include and exclude 
annihilation effects respectively. The inset details their behavior for small r. Notice that \(p(k r ,r)\ "heals" to (p(—iK,r) within 
0.1 fm away from the origin. 
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FIG. 11: Comparison of the bound state wave function obtained numerically (curve (a)) for the isoscalar one boson exchange 
potential of Eq. 11711 . with that given by the analytical solution, Eq. 1A.10II . for the equivalent Bargmann potential (curve (6)). 
Curve (c) shows the wave function in a gaussian potential adjusted to give a binding energy of 10 MeV |39l |. 
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FIG. 12: Variation of the spatial wave function ^>(0), Eq. I|A.10(I . with binding energy. 



